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1 Introduction and Some Preliminaries 

Selberg-type, positive definite symmetric matrices, gamma and beta integrals are considered by 
several authors; see Askey and Richards (1989), Gupta and Kabe (2005). Mathai (1997, p. 231, 
4.1.2) lists a positive signature symmetric matrix Selberg-type gamma integral. Recently Gupta 
and Kabe (2008) presented some results on the Selberg-type squared matrices gamma and beta 
integrals. They also covered Selberg-type skew symmetric matrices integrals. In this note, we 
extend the existing results in the literature for real normed division algebras to cover real, complex, 
quaternion and octonion spaces simultaneously. We mainly focus on Selberg-type square matrices 
Kummer-gamma and Kummer-beta integrals. 

Kabe (1984) developed the hypercomplex multivariate analysis distribution theory to consider 
real, complex, quaternion and octonion spaces simultaneously. Afterward many researchers ex- 
tended the results in many disciplines. Among them, the works of Prof. Jose A. Diaz-Garcia and 
his colleague should be acknowledged. The reader is refereed to Diaz-Garcia and Gutierrez-Jaimez 
(2009a, 2009b, 2011). 

A vector-space is always a finite-dimensional module over the field of real numbers. An algebra 
J 7 is a vector space that is equipped by a bilinear map m : J- x T — > T termed multiplication and 
a non-zero element 1 € J 7 termed the unit such that mil, a) = m(a, 1) = 1. As usual abbreviate 
m(a, b) = ah as ab. Do not assume J- is associative. 

An algebra J 7 is a division algebra if given a, b G J 7 , then ab = implied a = or b = 0. Equivalently, 
J 7 is a division algebra if the operation of left and right multiplications by any non-zero element 
is invertible. A normed division algebra is an algebra T that is also a normed vector space with 
\\ab\\ = [|o||||6||. This implies that J 7 is a division algebra and ||1|| = 1. 
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There are exactly four normed division algebras (Baez, 2002): 

(1) Real Numbers (E) 

(2) Complex Numbers (C) 

(3) Quaternions 

(4) Octonions 

moreover, they are the only alternative division algebras, and all division algebras have a real 
dimension of 1, 2, ,4 or 8, which is denoted by /3, see Baez (2002), Theorems 1, 2 and 3. The 
parameter a = 2/(3 is used, in other mathematical fields, see Edelman and Rao (2005). 
Let Cp : n be the linear space of all n x p matrices of rank n < p over T with m distinct positive 
singular values, where T denotes a real finite-dimensional normed division algebra. Let T ny -P be 
the set of all n x p matrices over T . The dimension of J rnx P over E is np(3. Let A G j^ ix p^ then 
A* = A T denotes the usual conjugate transpose. 

The set of matrices H\ G T ny -P such that i3"*ii"i = I p is a manifold denoted Vp, n , is termed 
the Stiefel manifold (Hi is also known as semi-orthogonal (f3 = 1), semi-unitary ((3 = 2), semi- 
symplectic (f3 = 4) and semi-exceptional type ((3 = 8) matrices, see Dray and Manogue, 1990). 
The dimension of V PjV , over E is (np(3 — \p(p + l)/3 — pj , is the maximal compact subgroup lA^(p) 
of Cp :P and consists of all matrices H £ J^p x p such that H*H = I p . Therefore, U^(p) is the real 
orthogonal group 0(p) ((3 = 1), the unitary group U(p) (f3 = 2), compact symplectic group S p (p) 
((3 = 4) or exceptional type matrices (p) ((3 = 8), for T = E, C,Q or O, respectively. 
Denote by Cp the real vector space of all S € J 7pxp such that S = S* . Let Bp be the cone of 
positive definite matrices S G jrpxp- then Bp is an open subset of Cp. Over E, Cp consists of 
symmetric matrices; over C, Hermitian matrices; over Q, quaternionic Hermitian matrices (also 
termed self-dual matrices) and over O, octonionic Hermitian matrices. Generically, the elements 
of are termed Hermitian matrices, irrespective of the nature of T . The dimension of Cp over 
E is (\\p{p + l)/3 + 2]). Let Vp 3 be the diagonal subgroup of C^ p consisting of all D G J rpxp , 
D = diag(di, . . . ,d p ). 

The surface area or volume of the Stiefel manifold Vp, n is given by 



Vol(V£ n ) = / {HldH^ 



np8 
2P1T2 



'Hievj!, n T p p (jfj ' 

and therefore 



p9 ( n£ 

dH l = }——(HldH{) = P - V (HUHi 

Vol(v£n) 1 2^ 1 1 

is the normalized invariant measure on Vp, n and dH, i.e. with (n = p), it defines the normalized 
Haar measure on U^(p) and Tp(a) denotes the multivariate gamma function for the space Cp, 
defined by 

T^dt \ f M|«-<tM-l , / A ^r A /r r / (i-l)/3\ 
T^(a) = / \A\ a 2 1 etr(-A)iiA = 7r 4 lir a-- —J, 
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where Re(a) > 1ZH_L2^ see Gross and Richards (1987). This can be obtained as a particular 
case of the generalized gamma function of weight for the space Cp with k = (k±, &2, . . . , k p ), 
k\ > &2 > ■ ■ • > k p > 0, taking k = (0, 0, . . . , 0) and for which Re(a) > (p — l)/3/2 — k p is defined 
by, see Gross and Richards (1987) and Diaz- Garcia and Gutiefrez-Jafmez (2009a) 

rJ(o,«)= / etr(-A)\Ar^-\ K (A)dA = {atrial 
where for A G Cp 

p-i 

q K (A) = \Ap\ k v\\\A^- k ^ 

1=1 

with A m = (a rs ), r,s = 1,2, ... ,m, m = 1, 2, ... ,p is termed the highest weight vector. Also (a)« 
denotes the generalized Pochhammer symbol of weight k, defined by 



(at = fl (a 



*, r?(«) i?(«) 

where Re(a) > (p — l)/3/2 — fc p and 
(a)j = a(a + 1) . . . (a + i - 1), 
is the standard Pochhammer symbol. 

The multivariate beta function for the space Cp is defined as (see Diaz- Garcia and Gutierrez- Jai'mez, 
2009a) 

B£(a,&) = / |X| a - 2 -il/p-XI 6 -^-^- 1 ^ 

Jo<x<i p 

= [ irr^-^ip + Yi-^dY 

JYeBp 

rg(°)rg(&) 
r^(a + &) ' 

where Y" = (I p — X)^ 1 - J p , Re(a),Re(6) > (p - l)/3/2. From Koev and Edelman (2006), 
(tr(X)) fc = ^Cf(X), 

where C«(X) denotes the zonal polynomials. 

Fix complex numbers ai, . . . , a r and b±, . . . , b s , and for all 1 < i < q and 1 < j < p do not allow 
~bi + (j — l)P/2 to be a non-negative integer. Then the hypergeometric function with one matrix 
argument r Fg is defined to be the real-analytic function on Cp given by the series 



■t a (ai,...,a r , 6 S , A) -2^2^77^ ~~k\ — ' 



{b x %...{b s f K 
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For convergence properties see Gross and Richards (1987). 
Then we have 

o#)=ee1t=^). 

k=0 k 

For more details and results the interested reader is referred to Diaz-Garcia and Gutierrez-Jafmez 
(2011), Dimitriu (2002), and Gross and Richards (1987). 



2 Matrix Variate Distributions 

In this study we consider some matrix variate distribution for our purpose. Consider the following 
definitions from Diaz-Garcia and Gutierrez- Jai'mez, (2009a, b) 

Definition 2.1. Let X G £p, n be a random matrix. 

1. (Matrix Variate Normal ) X is said to have a matrix variate normal distribution X ~ 

3 ( 

IXpV 



N^ x (/x, &), of mean /x and Cov(vec X*) = @ <g> £, if its density function is given by 



P 2 J P 

,„ , Ppn Bn Bp 

(2tt)~|5:|^"|0| — 



etrj-^E-^X-^r©- 1 ^-//)} 



2. (Wishart) Suppose that X has distribution X ~ N^ x (0, XI, I n ) and define S = X*X , then 
S is said to have a central Wishart distribution S ~ Wp (n, £) with n degrees of freedom and 
parameter S . Moreover, its density function is given by 



8pn 



Sl^^^etrj-^ n>(p-l)/3. 



ft 2 \a\ l3(n ~ p+1) 1 ± f P 

2^4 (^)|S|^? ' 



5. (T Type JJ) X zs said to /ia?/e a T Type / distribution X ~ T^ p (u, I p ), if its density function 
is given by 

r p ( 0(n±u) \ 

' + 2 , i/>p. 



vr^r^ (£) 

TTws distribution is also termed the T or Pearson Type VII matricvariate distribution. 

4- ( Gegenbauer Type II) X is said to have a Gegenbauer Type II distribution X ~ Gn,p{v, I p ), 
if its density function is given by 



I p -X*X\ iiJ ^ ±11 '\ i/>(p-l)/2. 



r /3 ( P{n+u) 
L P\ 2 

This distribution is known in statistical bibliography as the inverted T or Pearson Type II 
matricvariate distribution. 
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5. (T-Laguerre Type II Ensemble) X is said to have a T-Laguerre Type II ensemble distri- 
bution X ~ TLn, p {v), if its density function is given by 

' 8(n-p+l) , 0(n+v) 

|X|-Hs— l -i|j p + X| —, v>(p-l)/3, n>(p-l)p. 



R /3 ( 0p 0n\ 
D P \ 2 ' 2 J 

This distribution is also known as the Studentized Wishart distribution 

6. (Gegenbauer-Laguerre type II Ensemble) X is said to have a Gegenbauer-Laguerre type 
II ensemble distribution X ~ GIn.piy), if its density function is given by 

! dtn-p+l) -, /3(v-p + l) 

|x| £L_^_ 1|Jp _ x| _^_^ v > {p _ l) ^ n >(p-l)p. 



B 



.13 ( 0p 0n\ 
V \ 2 ' 2 ) 



The following definitions are analogous to the results of Nagar and Gupta (2002) and Nagar 
and Cardeno (2001) respectively for real normed division algebras. 

Definition 2.2. (Kummer-Beta) X is said to have a Kummer-beta distribution X ~ KB p (ai, a 2 , 5]), 
if its density function is given by 

K 1 (a 1 ,a 2 ,P,S)etr(-SX) \X\ ai a l \I p - X\ a2 a \ < X < I p 

where using Theorem 4-5 of Diaz- Garcia and Gutierrez- J aimez, (2009a), the normalizing constant 
is given by 



{^(ai^AS)}- 1 = / etr(-Sr)|Fr- i£ ^ M - 1 t/ p -y 

Jo<Y<I p 

= Bp 3 (ai,a 2 ) iFf {a\,a\ + a 2 ; -£). 



(p-1)/3 -i 



Definition 2.3. (Kummer- Gamma) X is said to have a Kummer-gamma distribution X 
KGp(at\,a 2 ,H), if its density function is given by 



K 2 (a 1 ,a 2 ,(3,V)etv(-VX)\X\ a i-^ }l - 1 \I p + X\- a \ X > 0, 



where using Theorem 3.9 of Diaz- Garcia and Gutierrez- J aimez, (2009a), the normalizing constant 
is given by 

{K 2 (a u a 2 ,p^)}~ 1 = [ etr(-Sr)|F| ai - i£ ^ 1£ - 1 |/ p + rr a2 dr 

Jy>o 

oo „ 

= EErf / C^-^Y)\Yr-^^\I p + Y\-^Y 
rTn^K' j Y> o 



k=0 k 



1 r^(ai,K)r^(a 2 -ai,-K) c p^_j^ 
"~" ^ ! if (a 2 ) 



fc=0 « 
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The following result is an extension to Theorem 3.1 of Nagar and Gupta (2002) for real normed 
division algebras. 

Lemma 2.1. Let U ~ K B$(ai, a 2 , £). Then for given matrices Q € Bp , (fi > the random 
matrix X defined by 

X = (ft- ¥)5E7(n - ¥)§ + * 

has generalized matrix variate Kummer-beta denoted by X ~ GKBp (ai, 02, 0, fi, ^) wi/i f/ie 
following density function 

Ci(ai,a 2 ,/3,e,fi,*)etr(-0X)|X - ^"i-*^" 1 ^ - Xl 02- * 2 ^ -1 , 

lu/iere *<X<fi ; = (fi- *)-s£(n - and 

Ci(ai, a 2 , P, 0, *) = Ki(ai, a 2 , 0, (O - *)^0(r2 - etr(0*)|O - 

Proof: The proof directly follows by Definition 3.2 and the fact that the Jacobian of transformation 

is j(u -»• x) = [n - ^i -1 ^- 1 . 

The following result is an extension to Theorem 2.2 of Nagar and Cardeno (2001) for real 
normed division algebras. 

Lemma 2.2. Let U ~ KGp(ai,a 2 ,Ti). Then for given matrices \I/,fi € Bp, f/ie random matrix 
X defined by 

X = (ft + *)it/(n + *)2 + * 

/ias generalized matrix variate Kummer-gamma denoted by X ~ GKGp(pt\, a 2 , 0, ft, wi/i i/ie 
following density function 

C 2 (ai,a2,/3,0,n,*)etr(-0X)|X - ^r 1-12 ^ -1 ^ + X|" a2 , 
w/iere * < X, = (fi + *)~5E(n + *)~§ and 

C 2 (ai,a 2 ,/3,0,",*) = K 2 ( ai ,a 2 ,f3, (ft + *)^0(ft + *)l)etr(0*)|n + *|- Ql+Q2 
Proof: The proof directly follows by Definition 3.3 and the fact that the Jacobian of transformation 

is j(u -> x) = |n + *|- £ii f i1 - 1 . 

3 Selberg Type Square Matrices Integrals 

In this section we are interested in evaluating the integrals of the form 

r p 

J h(A) rj(A* - XjfdA, (1) 

i<j 

where A = Diag(Ai, . . . , X p ), Ai > . . . > X p > 0, for some functions h. In this regard, we need 
the following essential result due to Diaz-Garcia and Gutierrez- Jafmez, (2009a). Note that similar 
results can be also found in Diaz-Garcia and Gutierrez- Jamiez, (2009b) for general form of the 
distributions defined in Definition 3.1. However the purpose of this study is not to find the joint 
density of eigenvalues, is to evaluate integrals of the form ([I]). 
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Theorem 3.1. Let X € Bp be a random matrix with density function f(X). Then the joint density 
function of the eigenvalues X±, . . . , X p of X is 

lpi/3+e p , 
g(X 1 ,...,X P ) = -rr^\H^- X rf f(HAH*)dH, 



P \ 2 ) »<J 

where dH is the normalized Haar measure and 

0, p = i 

-P, /? = 2 
-2p, /3 = 4 
-4p, (3 = 8. 

Lemma 3.1. Let X € Bp be a random matrix with density function f(X). Then we have 

r ( f * \ tS(% 

I [ f(HAH*) I](A, - A^dA dH = 

Proof: Since j A g(Xi, . . . , A p )dA = 1, by making use of Theorem 4.1 and changing the order of 
integration, the result follows. 

In sequel we proceed by giving some examples of gamma and beta integrals for real normed 
division algebras. 

3.1 Examples 

In this part, we give some examples of Selberg-type integrals. 

Examplel. Based on Definition 3.1, and the fact that H*H = I p , we have the following 
results using Lemma 4.1 

(1:1) Suppose that X ~ W p (n,I p ), then (proposed by Mathai, 1997) 



8pn 

2~ 



r /3 (p/3 
L P \ 2 



/ A |A|^-etr{-§A}n(A J -A^dA_ ^ ^ 



(1:2) Suppose that X ~ t£ :P (i/,I p ), then (proposed by Gupta and Kabe, 2008) 



/ \I P + A 2 \~^~ TT(Aj - XjfdA = — ^ ) / 

Ja P ' fjf 3> ^PP( P -n)+e T P(^^ 

(1:3) Suppose that X ~ Gn,p{v,I p ), then 

I \Ip-^\^^\{ { X l -X ] fdK = - 1 

■'A • , • IT 2 



i* ¥ £ 



t<3 
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(1:4) Suppose that X ~ TLn, P {y), then (proposed by Gupta and Kabe, 2005) 

r P 

/ |7V| — ^ + A| ^H^- A i) dA 



, . . , « . . ^ 2 2' 2 



(1:5) Suppose that X ~ GL^ !P (u), then 



P ( m \ rP ( §E §2L 

,£(2i^±l)_l lr . ,£(fiz£±i) tt m x x«^a _ P V 2 / P V 2 ' 2 



/ iAi^-i/ p -Ai-^n^-^ dA 




Example2. 

(2:1) Suppose that X ~ KB P (a±, a 2 , £), then by Definition 3.2 we get 

p 

etr (SHAH*) \I P - A\ ai - iE ^ 1Il - 1 \A\ a2 - iE ^ M - 1 TT(Ai - A^dAdJf 

I H J A iKj 

= i 2flT iif(ai,ai + a 2 ;-E)etr(E). (2) 

It is easily seen that taking S = 0, a\ = 4p and a 2 = in ([5]) gives item (1:5) in the above. 
Note that the integral in ([5]) is not of required form, thus making use of (see Diaz-Garcia and 
Gutierrez-Jahnez, 2009a) 



oo 

/ etr (SHAH*) dH = 



CkCEHAH*) 



k=0 re 



oo 



^^ g K (E)g K (A) 

k=0 k KV p; 



we obtain 



k=0 k ^ P> JA i<j 

(*f)BS(a lt a 2 ) 



L Ff(ai,ai + a 2 ;-£)etr(£). (4) 
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(2:2) Suppose that X ~ KGp (a± , ai , E), then by Definition 3.3 we get 

/ / etr(S^A//*)|A| ai ~ I ^ li! - 1 |I p + Ar a2 J](A; - A^dAdi* 

{^(ai.aa.^S)}" 1 . (5) 



r /3 (pft 
L P [ 2 



Similar result to (4) can be directly followed for (5) using (3). 
(2:3) Suppose that X ~ GKBp(ai,a2, @, fi, *), then by Lemma 3.1 we have 

/ / etr (-&HAH*) \HAH* - * 

JH J A 

P 

|n - HAH*\ a *- iE ^ }IL - 1 - Xj fdAdH 

TS(4)BS(a lia2 ) 



(p-1)/3 ± 



iiT(ai.ai + «2; -S) etr(-e*)|fi - -\ (6) 

where E = (fi - *)20(fi - 
(2:4) Suppose that X ~ GKGp 3 («i, a 2 , 0, fi, then by Lemma 3.2 we have 



/ / etr (-6 

JH J A 



QHAH*) \HAH* - 

p 

|« + HAH*\- a * - XjfdAdH 

i<j 

r /3 (pi 



P \ 2 



{CMai.aa.&e,^*)}- 1 . (7) 
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